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ABSTRACT 

Context. Water is a primordial species in the emergence of life, and comets may have brought a large fraction to Earth to form 
the oceans. To understand the evolution of water from the first stages of star formation to the formation of planets and comets, the 
HDO/H 2 ratio is a powerful diagnostic. 

Aims. Our aim is to determine precisely the abundance distribution of HDO towards the low-mass protostar IRAS 16293-2422 and 
learn more about the water formation mechanisms by determining the HDO/H 2 abundance ratio. 

Methods. A spectral survey of the source IRAS 16293-2422 was carried out in the framework of the CHESS (Chemical HErschel 
Surveys of Star forming regions) Herschel Key program with the HIFI (Heterodyne Instrument for the Far-Infrared) instrument, 
allowing detection of numerous HDO lines. Other transitions have been observed previously with ground-based telescopes. The 
spherical Monte Carlo radiative transfer code RATRAN was used to reproduce the observed line profiles of HDO by assuming an 
abundance jump. To determine the H 2 abundance throughout the envelope, a similar study was made of the H^O observed lines, as 
the H 2 main isotope lines are contaminated by the outflows. 

Results. It is the first time that so many HDO and H' 8 transitions have been detected towards the same source with high spectral 
resolution. We derive an inner HDO abundance (T > 100 K) of about 1.7 x 10 7 and an outer HDO abundance (T < 100 K) of about 
8 x 10~". To reproduce the HDO absorption lines observed at 894 and 465 GHz, it is necessary to add an absorbing layer in front of 
the envelope. It may correspond to a water-rich layer created by the photodesorption of the ices at the edges of the molecular cloud. 
At a 3cr uncertainty, the HDO/H 2 ratio is 1.4-5.8% in the hot corino, whereas it is 0.2-2.2% in the outer envelope. It is estimated at 
~4.8% in the added absorbing layer. 

Conclusions. Although it is clearly higher than the cosmic D/H abundance, the HDO/H 2 ratio remains lower than the D/H ratio 
derived for other deuterated molecules observed in the same source. The similarity of the ratios derived in the hot corino and in the 
added absorbing layer suggests that water formed before the gravitational collapse of the protostar, contrary to formaldehyde and 
methanol, which formed later once the CO molecules had depleted on the grains. 
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1. Introduction 

Water is one of the most important molecules in the solar system 
and beyond, but an unresolved question remains of how water 
has evolved from cold prestellar cores to protoplanetary disks 
and consequently oceans for the Earth's specific, but probably 
not isolated, case. In addition to being a primordial ingredient 
in the emergence of life, this molecule plays an essential role in 
the process of star formation through the cooling of warm gas. 
It also controls the chemistry for many species, either in the gas 
phase or on the grain surfaces. In cold dense cores, the gas-phase 
water abundance is low, less than about 10~ 8 (see for example 



Bergin & Snell 


2002 


much higher va 


ues in 



Caselli et al. 2010 1, while it can reach 



ues in the outflowi ng gas of the low-mass pro 



tostars environments (~10 -10 ; Liseau et al. 1996 Lefloch 



* based on HerschelfHIFl observations. Herschel is an ES A space ob- 
servatory with science instruments provided by European-led Principal 
Investigator consortia and with important participation from NASA. 



et al. 2010, Kristensen et al. 2010) before the solar-type proto- 
planetary systems are formed. 

In standard gas-phase chemistry, H2O forms through ion- 
molecule reactions that leads to H30 + , which can dissociatively 
recombine to form H2O (e.g., Bates|19 86 l Rodgers & Charnley 



2002 ), or through the highly endothermic reaction O + H2 — > OH 
+ H, followed by the reaction of OH with H 2 (fWagner & Graff 
[19871 iHoUehbach & McKeel[T989l [Aflanson et aL||2004| >. The 
former process is typical of diffuse cloud conditions, whereas 
the latter only works in regions with high temperatures, such as 
shocks or hot cores. It has been realized for around 30 years that, 
in cold and dense regions, H2O is formed much more efficiently 
on the grains through a series of reactions involving O and H 



accreted from the gas (e.g., |Tielens & Hagen 


1982 


Jones & 


Williams 1984[ Mokrane et al.|2009 Dulieu et a 


.2010 


Cuppen 


et al.||2010|). Near protostars, the grain temperature rises above 



1996 Fraser et al.|2001| l that increases the H 2 gas-phase abun- 



dance in the inner part of the envelope called hot core/corino for 
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a high-mass/low-mass protostar (Meln ick et al.|2000{|Ceccarelli| 
|et aL|1999||2000al >. 

Deuterated water is likely to form with a similar chemistry to 
that of water. Theoretically, the HDO/H2O ratio should be high if 
water has been formed at low temperature, i.e. on cold grain sur- 
faces, and low if it is a product of the photodissociation region or 
shock chemistry. Indeed, this is caused by the zero-point energy 



difference in the vibrational potential ( Solomon & Woolf|1973[ l. 
As the deuterated species have higher reduced mass than their 
undeuterated counterparts, the zero-point energy of the deuter- 
ated species is lower (the difference in energy between H2O 
and HDO is 886 K; |Hewitt et al.|2005) . Consequently, the en- 
richment of deuterated water with respect to its main isotopo- 
logue takes place at low temperature. The deuteration fraction 
observed in high-mass hot cores is typically HDO/H2O < 10" 3 



( |Jacq et al.|1990l |Gensheimer et al.|1996| |Helmich et alp 996 ), 
although higher values (~10 _/ ) have recently been found for 
Orion ( |Persson et al. [2007] Bergin et al.||2010 i. This ratio has 
also been estimated in the inner envelope of low-mass proto- 



stars such as IRAS 16293-2422 dParise et al.|2005b at about 3% 



NGC1333-IRAS4B (J0rgens en & van D ishoeck 201 0) with an 
upper limit of 0.06%, and NGC1333-IRAS2A ( |Liu et al.|2011) 
with a lower limit of 1%. Deuterated water has also been sought 
in ices towards several protostars, and has allowed |Dartois et"aT] 
(2003J and Parise et al. ( 2003) to obtain upper limits of solid 
] JDO/H 2 from 0.5% to 2%. 

It is important to know the HDO/H2O ratio throughout the 
protostar envelope in order to determine how this ratio is pre- 
served after the dispersion of the envelope, when a protoplan- 
etary disk is left over, from which asteroids, comets, and plan- 
ets may form. Obviously, this ratio is needed in order to evalu- 
ate the contribution of comets for transferring water in Earth's 
oceans (Morbidelli et al. 2000 Raymond et al.|2004 Villanueva 



et al. |2009) , directly associated with the emergence of life. It 
seems therefore crucial to determine how similar the observed 
HDO/H2O ratios in protostellar environments are to those ob- 
served in comets and in ocean water on Earth (~ 0.02 %, e.g., 
Bockelee-Morvan et al. 1998 Lecuyer et al. 1998 ). Recently, 



Hartogh et al. ( 201 1) have reported a D/H ratio in the Jupiter 



family comet 103P/Hartley2, originating in the Kuiper belt, of 
0.016, suggesting that some of Earth's water comes from the 
same comet family. 

IRAS 16293-2422 (hereafter IRAS 16293) is a solar-type pro- 
tostar situated in the LDN 1689N cloud in Ophiuchus at a dis- 
tance of 120 pc ( |Knude & Hogl[T998l |Loinard et aLl[2008) . It 
is constituted of two cores IRAS16293A and IRAS 16293B sep- 
arated by ~5", and the source IRAS16293A itself could be a 
binary system (Wootten 198 9). Several outflows have also been 
detected in this source (|Castets et al.||2001[ |Stark et al.||2004[ 
Chandler et"aT1[2003| |Yeh et al.||2008) . This Class proto- 



star is the first source where a hot corino has been discovered 
jCeccarelli et al.||2000b| |Cazaux et aLl[2003j |Bottinelli et~aT 
2004). It is also a well-studied case thanks to its high deuterium 
fractionation. For example, the abundance of doubly deuterated 
formaldehyde has been estimated between 3 and 16% of the 
main isotopologue abundance (|Ceccarelli et al. 1998 2001 



Methanol also shows a high deuterium fractionation: about 30 
+ 20% for CH 2 DOH, 6 + 5% for CHD 2 OH, and -1.4% for 
CD3OH ( |Parise et al.|2004) . More recently, |Demyk et al.| ( [20T0) 
have determined a methyl formate deuterium fractionation of 
~15%, and Bacma mi et al.| ( |2010| ) have concluded that there 
is a ND/NH ratio between 30% and 100%. Singly deuterated 
water in IRAS 16293 has been studied with ground-based tele- 
scopes by |Stark et aT1 ( |2004) and |Parise etaLlpOOl) . The former 



find a constant abundance of 3 x 10~ 10 throughout the envelope 
with the JCMT observation of the HDO lo,i-Oo,o fundamental 
line at 465 GHz alone, whereas the latter obtain an inner abun- 
dance (where T > 100 K) X m — 1 x 10~ 7 and an outer abun- 
dance X oM < 1 x 10~ 9 using four transitions observed with the 
IRAM^}30m telescope, as well as a JCMl^] observation at 465 
GHz. Using the water abundances determined from ISO/LW^] 
observations, which suffer from both a high beam-dilution and 
very low spectral resolution ( |Ceccarelli et al.|2000a) , |Parise et al.| 
( 12005) ) estimated a deuteration ratio HDO/H 2 of about 3% 
in the hot corino and lower than 0.2% in the outer envelope. 
Thanks to the spectral survey carried out with Herschel/HIFI to- 
wards IRAS 16293 in th e framework of the CHESS key program 
(Ceccareili et al. 2010), numerous HDO and H^O transitions 
have been observed at high spectral resolution, allowing an ac- 
curate determination of the HDO/H2O ratio in this source. Heavy 
water (D2O) has also been detected in IRAS 16293 with the 
Oo,o fundamental ortho transition at 607 GHz with Herschel/HIFI 
( Vastel et al.|2010) and the 1 1,0- 1 0,1 fundamental para transition 
at 31 7 GHz withJCMT JButner et al.|2007| ). From both transi- 
tions, Vastel et al. (2010) estimated a D2O abundance of about 2 
X 10" 11 in the colder envelope. 

The main goal of this paper is to determine the abundance of 
HDO throughout the protostar envelope, using new HDO col- 
lisional coefficients with H2 computed by Faure et al. ( 201 1) 



and |Wiesenfeld et aL (201 1 1, and combining the Herschel/HIFI, 
JCMT and IRAM data in this source. In Sect. 2, we present the 
observations and in Sect. 3, the modeling of the HDO and H^O 
emission distribution in the source using the radiative transfer 
code RATRAN ( |Hogerheijde & van der Tak]|2000") . In Sect. 4, 
we discuss the derived water deuterium fractionation, and con- 
clude in Sect. 5. 



2. Observations 

In total, thirteen HDO transitions have been detected towards the 
solar-type protostar IRAS 16293, nine with the Herschel HIFI in- 
strument, three with the IRAM-30m, and one with the JCMT. 
Three other transitions observed with HIFI, although not de- 
tected, have been used to derive upper limits. In addition, fifteen 
transitions of H^O have been observed with HIFI. The ortho- 
Hl s O l^o-lo.i fundamental line is clearly detected, four other 
transitions are tentatively detected, and ten others were not de- 
tected but their upper limits give constraints on the water abun- 
dance. Numerous H^O lines have been detected but not used 
since contaminated by the outflows. The ortho-H^O 1 i,o - lo,i 
fundamental line has also been detected. 



2.1. HIFI data 

In the framework of the guaranteed time Key Program CHESS 
( |Ceccarelli et al.||2010) , we observed the low-mass protostar 
IRAS 16293 with the HIFI instrument ( |de Graauw et al |[20T0) 
onboard the Herschel Space Observatory (Pilbr att et al.||2010| ). 
A full spectral coverage was performed in the frequency ranges 
[480 - 1142 GHz] (HIFI bands 1 to 5), [1481 - 1510 GHz] 
(HIFI band 6a), and [1573 - 1798 GHz] (HIFI bands 6b and 7a). 
The observations were obtained in March 2010 and February 
201 1, using the HIFI Spectral Scan Double Beam Switch (DBS) 



Institut de RadioAstronomie Millimetrique 
James Clerk Maxwell Telescope 

Infrared Space Observatory/Long Wavelength Spectrometer 
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fast-chop mode with optimization of the continuum. Twelve 
HDO lines and fifteen Hl s O lines have been observed in these 
bands. Table [TJlists, for all transitions, the observed parameters. 
The HIFI Wide Band Spectrometer (WBS) was used, provid- 
ing a spectral resolution of 1.1 MHz (0.69 km s~' at 490 GHz 
and 0. 18kms 1 at 1800 GHz) over an instantaneous bandwidth 
of 4x1 GHz. The targeted coordinates were ar 2 ooo = 16 /l 32 m 
22*75, ^2000 = - 24° 28' 34.2", a position at equal distance of 
IRAS 16293 A and B, to measure the emission of both com- 
ponents. The DBS reference positions were situated approxi- 
mately 3' east and west of the source. The forward efficiency is 
about 0.96 at all frequencies. The main beam efficiencies used 
are shown in Table Q] and are the recommended values from 
|Roeifsema eTaT1 ( (20T2l i. 

The data were processed using the standard HIFI pipeline 
up to frequency and amplitude calibrations (level 2) with the 
ESA-supported package HIPE 5.1 ( |Ott][20T0] > for all the bands 
except band 3a, processed with the package HIPE 5.2. In the 
selected observing mode, all lines were observed at least four 
times (if they are on a receiver band edge), but generally eight 
times (four in LSB and four in USB) for each polarization. To 
produce the final spectra, all observations were exported to the 
GILDAS/CLASS^f software. Using this package, the H and V 
polarizations at the line frequencies were averaged, weighting 
them by the observed noise for each spectra. We verified, for all 
spectra, that no emission from other species was present in the 
image band. HIFI operates as a double sideband (DSB) receiver, 
and the gains for the upper and lower sidebands are not neces- 
sarily equal. From the in-orbit performances of the instrument 
( Roelfsema et al.|2012[ ), a sideband ratio of unity is assumed 
for the HDO transition seen in absorption against the continuum 
(band 3b). The observed continuum is therefore divided by two 
to obtain the SSB (single sideband) continuum. The relative cali- 
bration budget error of the HIFI instrument is presented in Table 
7 of Roelfsema et al. ( 2012| l. Considering all the upper limits and 
estimated errors, we assumed an overall calibration uncertainty 
of 15% for each line. 



2.2. IRAM data 

The IRAS 1 6293 Millimeter A nd Submillimeter Spectral Survey 
(TIMASSS; |Caux et al.|20"TT| ) was performed between 2004 and 
2007 at the IRAM-30m telescope in the frequency range 80-280 
GHz and at the JCMT-15m telescope in the frequency range 328- 
368 GHz. The observations were centered on the IRAS 16293 B 
source at a 2000 = 16 /! 32'" 22 s .6, 6 20 oo= - 24° 28' 33". In this 
survey, four HDO lines were detected, but only three of them 
(80.578, 225.897, and 241.561 GHz) have been used in this pa- 
per, as the fourth transition at 266.161 GHz lies in a part of the 



surve y where the calibration uncertainty is very high ( Caux et al 



2011 



The spectral resolution is 0.31 MHz (~1.2kms r ) at 81 
GHz~and 1 MHz (-1.3 kms" 1 ) at 226 GHz and 242 GHz. 

2.3. JCMTdata 

The HDO lo,i-Oo,o fund amental transition a t 464.924 GHz was 
previously observed by Stark et al. (2004) on the IRAS 16293 
A source (a 20 oo = 16* 32 ffl 22 s .85, £ 2 ooo= - 24° 28' 35.5") and 



by |Parise et aE1p005) on the IRAS 16293 B source. The JCMT 
beam at this frequency is about 11". Despite the different point- 
ings, the profiles and the intensities of the line obtained are sim- 
ilar. This line shows both a wide emission and a narrow deep 



absorption (see Fig. [3J. The observations are not optimized for 
an accurate continuum determination. However, the continuum 
level is a crucial parameter in the following modeling. We see 
thereafter that we can estimate it at about 1 K thanks to the mod- 
eling and the dust emissivity model we chose. It consequently 
means that the narrow self-absorption completely absorbs the 
continuum. 



HDO level diagram 
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Fig. 1. Energy level diagram of the HDO lines. In red, IRAM- 
30m observations; in green, JCMT observation; in blue, HIFI 
observations. The solid arrows show the detected transitions and 
the dashed arrows the undetected transitions. The frequencies 
are written in GHz. 



The antenna temperatures of the observations were converted 
to the r m b scale, using the values of the main beam (B e ff) an d for- 
ward efficiencies (F e $) given in Table[T]with the usual relation: 



= T* A x 



(1) 



Table[T]and the energy level diagram in Fig.[T]summarize the 
HDO, H^O and H^ 7 transitions observed towards IRAS 16293 
with their flux or upper limit. The 3<x upper limits on the inte- 
grated line intensity are derived following the relation 

3cr (K.km s- 1 ) = 3 x rms x V2 x dv x FWHM (2) 

with rms (root mean square) in K, dv, the channel width, in 
kms -1 and FWHM (full width at half maximum) in kms -1 . 
We assume FWHM = 5kms _I , which is the average emission 
linewidth. The FWHM given in Table[T]was determined with the 
CASSIS^Jsoftware by fitting the detected lines with a Gaussian. 

Using the available spectroscopic databases JPL ( |Pickett| 
et al.| 1998) 1 and CDMS ( jMiiller et al.|200T][2005| ), we carefully 
checked that none of our lines are contaminated by other species. 



3. Modeling and results 

3.1. Modeling 

The spherical Monte Carlo ID radiative transfer code RATRAN 
( |Hogerheijde & van der T ak 2000 ), which takes radiative pump- 



http://www.iram.fr/IRAMFR/GILDAS/ 



5 CASSIS (http://cassis.cesr.fr) has been developed by IRAP- 
UPS/CNRS. 
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Table 1. Parameters for the observed HDO, H' s O, H' 7 Q, and HD ls O lines (I) . 
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1655.8676 


2i,2-lo,i 


114 (4) 


5.46 x 10- 2 


HIFI 6b 


12.8 


0.96 


0.71 


226 


< 1.66 




a 


o-H 2 8 


1719.2502 


3o,3"2l,2 


196 (4) 


5.12 x 10~ 2 


HIFI 7a 


12.3 


0.96 


0.71 


277 


< 2.04 




a 


p-H''0 


1107.1669 


1 1,1 -Oo.o 


53 


1.81 x lO" 2 


HIFI 4b 


19.1 


0.96 


0.74 


47 


< 0.34 




a 


o-H 2 '0 


552.0210 


1 i,o- lo.i 


61 (4 » 


3.37 x 10" 3 


HIFI la 


38.4 


0.96 


0.75 


11 


0.58 


7.1 


a 


HD ls O 


883.1894 


1 l.i-Oo.o 


42 


7.96 x 10" 3 


HIFI 3b 


24.0 


0.96 


0.75 


22 


< 0.16 




a 



(1) The frequencies, the upper state energies (E up ), and the Einstein coefficients (Ay) of HDO and H 2 s O come from the spectroscopic catalog JPL 
I Pickett et al. 1998} . The ortho/para separation of H 2 s O was carried out in the CASSIS database. 
U) The rms is computed for a spectral resolution of 0.7 km s~' for the HDO transitions and 0.6 km s~ 
(3) References : a) this work; b) |Caux et aT1 < |201 1[ >; c) |Stark et aLl(2004f ; d) |Parise et aT1j2005f 
<4) The o-H 2 8 and o-H 2 7 upper energy levels quoted in the table assume that the fundamental 1 1 
<5) This value corresponds to the ratio between the beam efficiency and the forward efficiency. 



for the H 2 



O, H 2 7 0, and HD 18 transitions. 



-1 i level lies at 34 K and not at K. 



ing by continuum emission from dust into account, has been em- 
ployed to compute the intensity of the molecular lines and the 
dust continuum. An input model describing the molecular H2 
density, gas temperature, and velocity field profiles is required, 
in order to define the spherical region into many radial cells. The 
code applies a Monte Carlo method to iteratively converge on 
the mean radiation field, J v . Level populations of HDO, H^O, 
H^O, and HD ls O can then be calculated, once J v is determined 
for every cell. These level populations are required to map the 
emission distribution throughout the cloud. 

The source structure used in the modeling was determined 
by Crimier et al. (2010) with a radius extending from 22 AU to 
6100 AU (see Fig. |2J; however, the structure in the hot corino 
region (T > 100 K) may be uncertain since disks probably ex- 
ist in the inner part of the protostellar envelopes. Nevertheless, 
since the disk characteristics are unknown, we keep the struc- 
ture as it is. The radial velocity v r - y/2GM/r (where M is the 



stellar mass, G the gravitational constant and r the radius) is cal- 
culated for a stellar mass of 1 M G . For a higher mass (~ 2 M ), 
the line widths become too broad to reproduce the line profiles. 
Moreover, the value of 1 M Q agrees with the mass of the core A, 
while the mass of the core B is at most about 0.1 M Q ( |Ceccarefii| 
|et al.|2000a| |Bottinelli et aLl[2004l |Caux et al.|2"0Tl) . For a ra- 
dius greater than 1280 AU, the envelope is considered as static 
(the velocity is fixed at 0). This radius corresponds to a change 
in the slope of the density profile, marking the transition of the 
collapsing/static envelope ( Shu|1977} , and its value has been de- 
termined by Crimier et al. ( 2010[ ). To reproduce the widths of 
the absorption lines, the turbulence Doppler b-parameter (equal 
to 0.6 x FWHM) is fixed at 0.3 km s . If a lower (respectively 
higher) value is adopted, the modeled absorption lines become 
too narrow (respectively too broad) compared with the observa- 
tions. As source A is more massive than source B, we consid- 
ered that the structure is centered on core A. According to in- 
terferometric data of the HDO 3i 2-22,1 line at 226 GHz, deuter- 
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Fig. 2. H2 density (solid line) and gas temperature (dashed line) 
structure of IRAS 16293 determined by |Crimier et al.| ( |2"0T0) . 



The expected abundance profile of water is added (dashed-dotted 
line) on an arbitrary Y-scale. In the colder envelope (T<100K), 
the water molecules are trapped on the grain mantles, whereas in 
the inner part of the envelope, they are released in gas phase by 
thermal heating, leading to an enhancement of the abundance. 
In the outer part of the molecular cloud (Ay ~ 1-4), the water 
abundance can also increase by photodesorption mechanisms. A 
freeze-out timescale much longer than the protostellar lifetime 
can also lead to an enhanced abundance when the density is low 
(~ 10 4 ■ 



cm J ). 



ated water is mainly emitted by core A ( J0rgens en et al.|201 1 1. 
Therefore, the assumption of a ID modeling centered on source 
A seems quite reasonable. Because the IRAM observations were 
pointed on the IRAS 16293 B source and not on core A and the 
beam is quite small (~ 11") at 226 and 242 GHz, we carefully 
convolved the resulting map with the telescope beam profile cen- 
tered on core B to get the model spectra. 

To fit the continuum observed with HIFI from band 1 to band 
4, the dust opacity as function of the frequency has to be con- 
strained by a power-law emissivity model : 



K = K0\ 



\ v 0) 



(3) 



with /? = 1.8, kq = 15 cm 2 /gdust) an d v o = 10 12 Hz, which has 
been used as input in the RATRAN modeling. Using this emis- 
sivity model and the source structure described above, the con- 
tinuum of the HDO lo,i-Oo,o fundamental line at 465 GHz is pre- 
dicted at ~1 K. The continuum determination at this frequency 
is essential for the modeling because this line shows a deep 
and narrow self-absorption. The modeled continuum is shown 
in Figs. [3] and |4] for the two HDO absorption lines, whereas it 
has been subtracted for all the lines that only present emission. 
The profiles obtained are resampled at the spectral resolution 
of the observations. Smoothing was applied on some observa- 
tions when the line is undetected or weakly detected (see Figs. 
|4]and 10 1. For the determination of the best fit parameters by^ 2 
minimization (see Sect. 3.3), the spectra are also resampled to 
a same spectral resolution for all the lines (0.7 km s -1 for HDO 
and 0.6kms _1 for H^O, H^O, and HD ls O, which correspond 
to the lowest spectral resolution of the HIFI data). 

Previous studies (e.g. van K empen et al.|2008) have shown 
that in protostellar envelopes, the dust can become optically 



thick, preventing water emission from deep in the envelope to 
escape. The dust continuum could therefore hide the higher fre- 
quencies HDO and H^O lines coming from the hot corino, de- 
pending on the source structure and the dust emissivity model 
chosen. In fact, we notice here that the dust optical depth is al- 
ways lower than 1 for all the frequencies lower than 1500 GHz. 
For the transitions quoted in Table[T]with frequencies higher than 
1500 GHz, the dust opacity only exceeds 1 for temperatures well 
over 100 K. Consequently, in the case of IRAS 16293, the dust 
opacity should not hide the higher frequencies HDO and H^O 
lines as the dust opacity only becomes thick in the very deep part 
of the inner envelope. 



3.2. HDO collisional rate coefficients 

The HDO collisional rates used in this study have recently been 
computed by |Faure et al.| ( |2011| l for para-H2(/2 = 0, 2) and 
ortho-H2 (72 = 1) in the temperature range 5-300 K, and for 
all the transitions with an upper energy level less than 444 K. 
The methodology used by Faure et al. (2011) is described in 
detail in |Scribano et al.| ( |2010[ ) and [Wiesenfeld et al.| ( |2011) . 
These authors present detailed comparisons between the three 
water-hydrogen isotopologues H2O-H2, HDO-H2, and D2O- 
H2. Significant differences were observed in the cross sections 
and rates and were attributed to symmetry, kinematics, and in- 
tramolecular (monomer) geometry effects. Moreover, in the case 
of HDO, rate coefficients with H2 were found to be significantly 
larger, by up to three orders of magnitude, than the (scaled) H2O- 
He rate coefficients of Green (1989), which are currently em- 
ployed in astronomical models (see Fig. 2 of |Faure et al.|20lT) . 
A significant impact of the new HDO rate coefficients is thus 
expected in the determination of interstellar HDO abundances, 
as examined in Sect. 3.3. In the following, we assume that the 
ortho-to-para ratio of H2 is at local thermodynamic equilibrium 
(LTE) in each cell of the envelope. The collisional rates for para- 
H2 were also summed and averaged by assuming a thermal dis- 
tribution of J2 = 0, 2. 

3.3. Determination of the HDO abundance 

In a first step, we ran a grid of models with one abundance 
jump at 100 K, in the so-called hot corino ( |Ceccarelli et al.|1996 



2000b). The abundance of water is expected to be higher in the 
hot corino than in the colder envelope, since the water molecules 
contained in the icy grain mantles are released in gas phase, 
when the temperature is higher than the sublimation tempera- 
ture of water, ~100K (Frase Fet al.|200l) . Thereafter, the inner 
abundance (T > 100K) will be designated by X\ n and the outer 
abundance (T < 100 K) by X out . Both are free parameters and 
their best fit values are determined by a x 2 minimization. To 
take the line profile into account, the x 1 is computed from the 
observed and modeled spectra resampled at a same spectral res- 
olution (0.7 km s for HDO and 0.6 km s _1 for Hj 8 0) according 
to the following formalism: 



N nam 



X 



obsjj 



-r, 



mod.ij 



-f pf rmsf + (Cak X 7/ obs ,ij) 2 



(4) 



with N the number of lines i, n c h an the number of channels j for 
each line, r o b s ,;j and T mo d,ij the intensity observed and predicted 
by the model respectively in channel j of the line i, rms[ the 
rms at the 0.7kms~' (or 0.6kms _I ) spectral resolution (given 
in Table 1), and Cal\ the calibration uncertainty. We assumed an 
overall calibration uncertainty of 15% for each detected line. 
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The best fit result for this grid of models gives an inner abun- 
dance X m — 1.9 x 10~ 7 and an outer abundance X out = 5 x 
1CT 11 with a reduced^ 2 of 3.2. However this model (as well as 
the other models of this grid) predicts absorption lines that are 
weaker than observed (see Fig. [3]). To correctly reproduce the 
depth of the absorption in the 465 and 894 GHz lines, while not 
introducing extra emission in the remaining ones, it is necessary 
to add an absorbing layer in front of the IRAS 16293 envelope. 
We notice that the results mainly depend on the HDO column 
density of the layer, and they are insensitive to the density (for 
densities lower than ~ 10 s cirT 3 ) and to the temperature (for 
temperatures lower than ~ 30 K). The HDO column density of 
the layer must be about 2.3 x 10 13 cirT 2 . To understand from 
where the HDO absorptions arise, we can consider two different 
options. 

- First, the source structure should be extended to a higher 
outer radius to include the whole molecular cloud. If the ab- 
sorbing layer comes from the surrounding gas of the molec- 
ular cloud harboring IRAS 16293 and if we assume that the 
HDO abundance remains constant as in the outer envelope, 
this would imply a surrounding H2 column density of about 
2.9 x 10 23 cm -2 . But this value is considerably too high 
for a molecular cloud. The N(Ha) column density of t he p 



HDOI^-Oqo 893.63869 GHz 



HDO 1 01 -0 , 464.92452 GHz 



Oph cloud has been e stimated at ~ 1.5 X 10 22 cm 2 by 



Dishoeck et al. 



( |1995| | and at > 5 x 10 22 cirT 2 by 



Caux et al. 



( 1999). Adding a molecular cloud with these characteristics 
to the structure determined by |Crimier et al.|f2010| l does not 
allow deep absorption lines to be modeled. This hypothesis 
is therefore insufficient to explain the deep HDO absorption 
lines. 

- Consequently, the only way to explain the absorption lines 
consists in assuming a drop abundance structure. Such a 
structure in the low-mass protostars has already been in- 
ferred for several molecules like CO and H2CO (Schoier 
et al.|2004| |j0rgensen et al.|2004[ [2053) . The raising of the 



abundance in the outermost regions of the envelope is ex 
plained by a longer depletion timescale than the lifetime of 
the protostars (~10 4 -1 5 years) for H? densities lower than 



10 4 -10 5 cirT 3 (e.g., |Caselli et al.||l999[ |j0rgensen et al 



2004). On the other hand, |Hollenbach et al.| ( |2009| ) have" 
shown that, at the edge of molecular clouds, the icy man- 
tles are photodesorbed by the UV photons, giving rise to an 
extended layer with a higher water abundance for a visual ex- 
tinction Ay~ 1-4 mag (for Go = 1). Assuming an Ay of ~l-4 
mag and the relation N(H 2 )/A V = 9.4 x 10 20 cirT 2 mag -1 
( Frerking et al.|1982] l, the abundance of HDO in this water- 
rich layer created by the photodesorption of the ices should 
be about 6 x 10~ 9 - 2.4 x 10~ 8 . We see thereafter, when an- 
alyzing the Hl s O lines, that this hypothesis is nicely con- 



sistent with the theoretical predictions by Hollenbach et al. 
(2009). For water, both effects may play a role. 

As the absorption of the two fundamental lines is reproduced 
with this outer component (see Fig. [3}, we ran a second grid of 
models adding this absorbing layer. The best fit is obtained for 
X in = 1.7 x 10~ 7 and X out = 8 x lO" 11 with a x 2 ed of 2.4. With the 
current parameters, a simultaneous fit of all the HIFI, IRAM, and 
JCMT data is presented in Fig. [4] The contours of^ 2 at lcr, 2<r, 
and 3cr, which respectively represent a confidence of 68.26%, 
95.44%, and 99.73% of enclosing the true values of X m and 
X out , are shown in Fig. [5] computed with the method of |Lampton| 
et al. ( 1976 1. The contours at lcr, 2cr, and 3<x correspond respec- 
tively to~p=y 2 • +2.3, y 2 =y 2 ■ +6.17, and y 2 =y 2 . +11.8 when 

J A A mm ' A A mln ' A A mln 

the number of adjustable parameters is two (here X m and X out ). 





-10 -5 5 10 15 

v LSH (km s ') 



Fig. 3. In black: HDO 1 i,i-Oo,o and li,n-Oo.o absorption lines ob- 
served at 894 GHz with HIFI and at 465 GHz with JCMT, re- 
spectively. In red: HDO modeling without adding the absorbing 
layer. In green: HDO modeling when adding an absorbing layer 
with a HDO column density of ~ 2.3 x 10 13 cm~ 2 . The con- 
tinuum for both the 894 GHz and 465 GHz lines refers to SSB 
data. 
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Fig. 5. x 2 contours at lcr, 2cr, and 3cr obtained when adding 
an absorbing layer with a HDO column density of 2.3 X 10 13 
cirT 2 to the structure. The results are obtained with the colli- 
sional coefficients with ortho and para-H2 determined by Faure 
|et al. ( |201 1) . The best-fit model is represented by the symbol 
"+". The reduced^ 2 is about 2.4. 



The confidence intervals of X m and X out at 3cr are 1.53 x 10 7 - 
2.21 x 10~ 7 and 4.6 x lO" 11 - 1.0 x 10~ 10 , respectively. 

In the hot corino, most of the HDO lines are optically thick, 
except the 8 1 , 226, 242, and 482 GHz transitions that show an 
opacity < 1. Since the emission from most of the HDO lines 
originates both in the hot corino and the outer envelope, the use 
of a simple rotational-diagram analysis to estimate the total col- 
umn density of HDO is not appropriate. Figure |6]represents the 
normalized populations of some levels (up to 3 1,2) as a function 
of the distance from the center of the protostar, computed by 
RATRAN through the equation of statistical equilibrium. These 
level populations are then used to determine the emission distri- 
bution of deuterated water by using a ray-tracing method. This 
figure also shows that the absorbing layer is mainly constrained 
by the ground state transitions as l^i-O^o, and lo,i-On,o- 

The 81 GHz 1 1 ,0 - 1 1,1 transition observed at IRAM-30m is 
not reproduced well by the best-fit model. The disagreement 
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Fig. 4. In black: HDO lines observed with HIFI, IRAM, and JCMT. In red: best-fit model obtained when adding an absorbing layer 
with an HDO column density of ~ 2.3 X 10 13 crrT 2 to the structure (see details in text). The best-fit inner abundance is 1.7 x 10~ 7 
and the best-fit outer abundance is 8 x 10 -11 . The continuum shown for both the 894 GHz and 465 GHz lines, refers to SSB data. 



could be due to a blending effect at this frequency. Nevertheless, 
it is ruled out for the species included in the JPL and CDMS 
databases. Calibration could also be questioned at this frequency. 
But it seems rejected since the profile obtained with the best- 
fit model differs from the observing profile. Another explana- 
tion could come from the assumption of spherical symmetry. 
Interferometric data of the HDO 2>ix^2,i transition show that, 
in addition to the main emission coming from the source A, a 
weaker component of HDO is emitted at about 6" from core A 
(see Fig. 20 in J0rgense n et al.| ( |2011| >). The beam size of the 
IRAM-30m telescope at 81 GHz encompasses this component 
and could explain the lack of predicted flux as well as the differ- 
ent line profile. This hypothesis seems consistent with the lower 
frequency HIFI lines (491, 509, and 600 GHz) for which the 
best-fit model misses flux. 

The derived abundances for HDO agrees with the 3cr upper 
limit (0.16Kkms _1 ) of the HD ls O 1 i,i-0 ,o transition observed 
in o ur spectra at 883 GHz. A ssuming an HD 16 0/HD 18 ratio of 
500 ( |Wilson & Rood|1994| l, value at the galactocentric distance 
of IRAS 16293, the modeling predicts an integrated intensity of 
0.15 Kkms- 1 . 



To obtain consistent profiles of the model compared with the 
observations, the LSR velocity, Vlsr, is about 4.2kms~' for the 
HIFI and JCMT lines, whereas it is about 3.6 km s _1 for the three 
IRAM lines. This discrepancy of the velocity in the modeling 
could be explained by the origin of the emission of the lines. 
The 4.2kms~' velocity component could represent the velocity 
of the envelope. For example, the velocity of the absorption lines 
of the fundamental transitions of D2O tracing the cold envelope 
of IRAS 16293 and detecte d by JCMT ( [Burner et al.|[2007] > and 
by HIFI ( |Vastel et al.|2010) > is 4.15 and 4.33 km s" 1 , respectively. 
Even if some HIFI lines trace both the inner part and the outer 
part of the envelope, the velocity of the outer envelope should 
dominate because of a smaller line widening in the outer enve- 
lope. In contrast, the flux of IRAM lines at 226 and 242 GHz 
only trace the hot corino (see Figs. [T] and [SJl. Their velocity (~ 
3.6kms~') is lower, in agreement with the velocity of the core 
A (y LS R -3.9 km s- 1 ). 

As previous studies regarding HDO used the collisional co- 
efficients computed by Green ( 1989 ), we have run a grid of mod- 
els with these rates for a comparison. The simultaneous best- 
fit of all the transitions using these rates is obtained for an in- 
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Fig. 6. Normalized level populations of HDO computed by 
RATRAN, as a function of the radius of the protostar envelope. 
To avoid confusion, the population of levels 2o,2, and 3i2 
are indicated by dashed lines. The HDO level population 
jump visible at ~1300 AU is an artifact due to the change of the 
velocity field at the interface between the infalling envelope and 
the static envelope. 



ner abundance X m = 2.0 x 10 and an outer abundance X out 
= 1 x 10 -11 . The inner abundance is similar to what we found 



above with the collisional coefficients determined by Faure et al. 
(201 1 ). However, the outer abundance is a factor of 8 lower. The 
X 1 also gives a higher value (xLa - 2-8), therefore showing that 
the observations are better reproduced with rates with H2 than 
with He. 



3.4. Determination of the water abundance 

To determine the H2O abundance throughout the envelope, we 
used all the H^O transitions in the HIFI range (see Table [TJi. 
The only transition observable in the TIMASSS spectral survey 
is contaminated by the CH3OCH3 species at 203.4 GHz. The 
profiles of the observed H^O lines suggest, by the presence of 
wings, that the lines are contaminated by the outflows. Unlike 



the H^O lines that show a width > lOkms 1 



the H' s O line 

widths are similar to the HDO lines (see Fig. [7] ~ 5 kms' 1 ) and 



no wings are seen in the H 2 O detected transitions (see Fig. 10 1. 
Consequently, the outflow does not probably contribute to a large 
extent to the emission of the H^O lines. 

We used the H?Q collisional coefficients determined by 
|Faure et al.|(|2007|l and assumed a H3, 6 0/H^ 8 ratio equal to 500 
(Wilson & Roodl 1994b, as well as a H^O ortho/para ratio of 3. 
As tor HDD, in a first step, we ran a grid of models for different 
inner and outer abundances without adding the absorbing layer. 
The best-fit parameters of the x 2 minimization are Xjj,(Hl O) = 
1 x 10~ 8 and X olU (H' 8 0) = 3 x 10" 11 and the reduced x 2 is about 
1.9. The inner abundance of water is therefore 5 x 10~ 6 , whereas 
the outer abundance is 1.5 x 10~ 8 . Figure [8] shows the^ 2 con- 
tours at lcr, 2<r, and 3cr. A calibration uncertainty of 15% has 
been assumed for the 548 GHz transition and the tentatively de- 
tected ones (1096, 1102, 994, and 745 GHz). At 3cr, the outer 
abundance of H^O varies between 9 x 10~ 12 and 4.9 x 10~ n , 
whereas the inner abundance is in the interval 7.6 x 10~ 9 - 2.1 x 
10~ 8 . Using the best-fit values found for X m and X out , it is nec- 




Vlsr (km s 1 ) 



Fig. 7. Comparison of the profiles of the HDO (in black) and 
Hj 8 (in red) 1 i,o-lo,i transitions at 509 and 548 GHz, respec- 
tively. 
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Fig. 8. x 1 contours at lcr, 2cr, and 3cr obtained for Hi 8 0. The 
best-fit model is represented by the symbol "+". The reduced ;f 2 
is about 1.9. 



essary, in order to reproduce the weak absorption of the 1 i,o-lo,i 
line at 548 GHz, to add an absorbing layer with a H^O column 
density of about 1 x 10 12 crrT 2 . For higher column densities, the 
models predict absorptions that are too deep for the ortho and 
para fundamental transitions (see Fig.|9]l. For a visual extinction 
Ay~ 1-4 mag, the H2O abundance in the photodesorption layer 
is 1.3 - 5.3 x 10~ 7 . This value is in very good agreement with 
the values predicted in photodesorption layers by the model of 
Hollenbach et al.|(|2009|, about 1.5 - 3 x 10~ 7 . The model predic- 

We cannot 



10 



tions with this absorbing layer are shown in Fig. 
rule out a remnant contribution of the outflows in view of the 
profile of 548 GHz line, but this outflow contribution to the bulk 
of the emission is very likely negligible for the other transitions. 

In the hot corino, among the o-Hi 8 O transitions, only the 489 
and 1181 GHz lines are optically thin, whereas among the p- 
H' 8 transitions, only the 970, 1189, and 1606 GHz lines show 
an opacity lower than 1. The emission of most of the transi- 
tions comes from both the hot corino and the outer envelope (see 
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Fig. 10. In black: H l 2 O lines observed with HIFI. In red: best-fit model obtained when adding an absorbing layer with a H\ 
column density of 1 x 10 12 cm -2 . The inner abundance is 1 x 10~ 8 and the outer abundance is 3 x 10 -11 . The detected line at a 
velocity of 15kms 1 in the panel of the 489 GHz transition is CH3OH at 489.0368 GHz. Also, the line observed at a velocity of 
-lOkms- 1 in the panel of the 970 GHz transition is p-H 2 Q at 970.3152 GHz. 



Fig. [9}. Similarly to HDO, a rotational diagram analysis there- 



fore cannot be considered for Hi 8 0. 



Hj O 1 i,o-lo, 1 transition at 552 GHz has also been 



The ortho 

detected. To check the validity of the results, we ran a model sim 
ilar to what has been done for Hj 8 0, assuming a ratio 18 0/ 17 of 
4 ( Wouterloot et al.|2008 1. Figure 11 shows the predicted mod- 
els for the 552 GHz line, as well as the para-H^O lij-Oo.o fun- 
damental line undetected at 1107 GHz for the best-fit of H' s O, 



i.e. for an inner 



H^O abundance of 2.5 



10 9 and an outer 



H^O abundance of 7.5 x 10~ 12 . A prediction with twice the in- 
ner abundance, situated in the 3<x contour of H^O, is also pre- 
l. We see that the H^O predictions agrees with 
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sented (Fig. 

the H^ 8 results, confirming the water abundances derived. 



4. Discussion 

The study of HDO in the deeply embedded low-mass protostar 
IRAS 16 293 had already bee n undertaken by |Stark et aT] ( [2004)i 
then by |Parise et al.| ( |2005| ), leading to different results, j Stark 



etal.| ( |2004t used a constant [HDO]/[H 2 ] abundance of 3 x 10" 10 
throughout the envelope to fit the lo,i-Oo,o transition at 465 GHz 
best, resulting in an HDO/H2O ratio in the warm inner envelope 
of a few times 10~ 4 . However, their best fit cannot reproduce the 
deep absorption feature. Also, unlike our study, they do not suc- 
ceed in reproducing the line width without introducing an out- 
flow contribution; however, the model obtained with the outflow 
component fails to match simultaneously the line width and the 
peak intensity. This could be explained by their source structure 
and the assumption of a constant abundance throughout the en- 
velope. The dust temperature is restricted to a range from 12 K 
to 1 15 K. A structure with a lower inner radius (consequently an 
higher temperature) and an abundance rising in the inner part of 
the envelope may allow the line to become broader and repro- 
duce the observation. Moreover, the peak intensities predicted 
by their model for the lines at 226 and 241 GHz (0.001 K and 
0.006 K, respectively) are clearly lower than what we observe 
here (0.30 and 0.32 K, respectively). Afterwards, the results ob- 
tained by Stark et al. ( 2004) have been questioned by the study of 



Parise et al.| ( |2005| l. Using IRAM and JCMT observations, they 
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1000 
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Fig. 9. Normalized level populations of Hi O computed by 
RATRAN, as a function of the radius of the protostar envelope. 
To avoid confusion, the population of levels lio and 3o,3 are in- 
dicated by dashed lines. 



H 2 17 1,0-1,,, 552.021 GHz 



H 2 ,7 1,,-Onc 1107.167 GHz 





Fig. 11. In black: Fundamental Hi 7 lines observed with HIFI. 
In red: model obtained with an inner H^O abundance of 2.5 x 
1(T 9 , an outer H' 7 abundance of 7.5 x 1(T 12 , and an H' 7 
column density in the absorbing layer of 2.5 X 10 11 cm 2 . In 
green: model obtained with an inner H l 2 7 abundance of 5.0 x 
!Q-\ an outer H' 7 abundance of 7.5 x 10" 12 , and an H' 7 Q 



-2 



column density in the absorbing layer of 2.5 X 10 11 cm 



have found an enhancement of the abundance of HDO in the hot 
corino (X[ n ~ 1 x 10~ 7 ) with respect to the outer envelope (X out 
< 1 x 10~ 9 ), in agreement with our results using both the colli - 
sional coefficients with He by Green ([T989 ) and with ortho and 
para-H2 by |Faure et al.| ( |20lT] ) (see Sect. 3.3). However, their re- 
sults did not take the line profiles into account, specifically the 
deep absorption of the ground state transitions, and only allowed 
a rather high upper limit of the outer abundance to be derived. 



Parise et al. (2005} conclude there is a jump in the HDO/H2O 



ratio in the inner part of the envelope, using the H2O abundance 
determined by Ceccarelli et al. ( 2000a| l. Nevertheless, our re- 
sults on the water abundance are quite different. Using H^O 
data from ISO/LWS, highly diluted in a beam width of about 
80", |Ceccarelli et al.| (|2000a) obtained an outer abundance of 
5 x 10"', while we found an outer abundance lower than 3.5 
X 10~ 8 . However, their derived inner abundance of ~ 3 x 10~ 6 
is not so very different from our value, X; n (H20) ~ 5 x 10~ 6 . 
The disagreement on determining the outer abundance can eas- 
ily be explained by the fact that |Ceccarelli et al.| ( |200 0a) used 



H\ data, contaminated by the outflows that must be added in 
the modeling. The gas temperature was computed by assuming 
the water abundance derived by the ISO observations. However, 
since the gas is practically thermally coupled with the dust, the 
impact in the gas temperature is certainly negligible. The largest 
difference between the gas and dust temperature is less than 10% 
( |Crimier et al.|2010) . The discrepancy on the water abundance 
therefore results in a different deuteration ratio in the outer en- 
velope. Considering the 3cr uncertainty, the enhancement of the 
water deuterium fractionation in the hot corino cannot be con- 



firmed in our study, in contrast to what was concluded in |Parise 
|et al.| p005] l. Table [2] summarizes the HDO/H 2 ratio deter- 
mined by our analysis in the inner and outer envelope, as well 
as an estimation of this ratio in the absorbing layer. With a 3cr 
uncertainty, the water fractionation could be similar throughout 
the cloud with a HDO/H 2 ratio between 1.4% and 5.8% in the 
inner part and between 0.2% and 2.2% in the outer part. Also 
the HDO/H2O ratio obtained here is in agreement with the up- 
per limits of solid HDO/H 2 (from 0.5% to 2%) determined by 
observations of OH and OD stretch bands in four low-mass pro- 
tostars ( |Dartois et al.|2003t|Parise et al.|2003) . 

A consequential result of this paper is the similarity of the 
HDO/H 2 ratio derived in the hot corino (~ 1.4-5.8%) and in the 
added outer absorbing layer (~ 4.8%). Indeed, the ratio is within 
the same order of magnitude, although the densities of these 
two regions are considerably different. The density in the hot 
corino is a few times 10 8 cm"" 3 . In the absorbing layer, it is about 
10 3 -10 5 cm" 3 . Consequently, water shows a different behav- 
ior from other molecules such as methanol and formaldehyde, 
that need CO ices to be formed. In particular, Bacmann et al. 
( |2003l |2007> show that the deuteration of H 2 CO and CH 3 OH 
increases with the CO depletion in starless dense cores, increas- 
ing itself with the H2 density (Bacman n et al.|2002 i. The H2CO 
and CH3OH deuteration is therefore sensitive to the density of 
the medium in which they form. In contrast, the HDO/H2O ra- 
tio does not show any difference for different H2 densities. This 
would therefore mean that water has formed before the collapse 
of the protostar and that the HDO/H2O ratio has been preserved 
during the gravitational collapse. The deuteration fractionation 
of water would therefore remain similar, both in the inner re- 
gion of the protostar, where the density has strongly increased, 
and in the outer region that has not been affected by the col- 
lapse. Another argument emphasizes this hypothesis. To obtain 
such high deuteration ratios of HDCO/H2CO (-15%; Loinard 
|eTaT1[2WT ) and CH 2 DOH/CH 3 OH (-30%; |Parise et al.||2004] i 
in IRAS 16293 compared with HDO/H 2 0~3%, the density at 
which the molecules form should be higher for the H2CO and 
CH3OH formation than for water formation. The collapse should 
already have started to allow the CO molecules to freeze-out and 
form formaldehyde and methanol at the grain surface. On the 
contrary, water would form at low densities in the early stages 
of the star formation before the protostellar collapse, as sug- 
gested by Dartois et al.| ( |2003) and Parise et al. ( 2003) >. This is 



also consistent with the fact that H2O ices appear at relatively 
low extinction in the direction of the Taurus cloud (e.g. Jones & 



Williams 1984). Recently, similar conclusions have been men- 



tioned in (Cazaux et al. (2011). Using a grain surface chemistry 
model, they show that the deuteration of formaldehyde is sen- 
sitive to the gas D/H ratio as the cloud undergoes gravitational 
collapse, while the HDO/H2O ratio is constant as the cloud col- 
lapses and is set during the formation of ices in the translucent 
cloud. 

The HDO/H 2 ratio in the low-mass protostar IRAS 16293 
is close to what is found in the low-mass protostar NGC1333- 
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Table 2. HDO/H 2 ratio 





Hot corino 


Outer envelope 


Photodesorption layer 


Best fit 


3a- 


Best fit 


3a- 


A v ~ 1 - 4 mag 


HDO 


1.7 x 10~ 7 


1.5- 2.2 x 10"' 


8 x 10" u 


4.6- 10.0 x 10" 11 


~ 0.6 - 2.4 x 10" 8 


H 2 


5 x 10- 6 


3.8 - 10.5 x 10- 6 


1.5 x 1(T 8 


4.5 - 24.5 x lQ- g 


~ 1.3- 5.3 x 10~ 7 


HDO/HoO 


3.4% 


1.4% - 5.8% 


0.5% 


0.2% - 2.2% 


~ 4.8%" 



" This ratio remains valid if the absorbing layer is not due to the photodesorption, but is the result of a freeze-out timescale longer than the 
protostellar age at low densities. 



IRAS 2 A, as determined by |Lme t al. (2011): higher than 1% 
in the hot corino and between 0.9% and 18% at 3cr. However, 
this ratio of a few percent does not seem typical of all the Class 
protostars since an upper limit of 6 x 10~ 4 has been deter- 
mined in the inner part of the envelope of the low-mass proto- 
star NGC1333-IRAS4B (J0r gensen & van Dishoeck|2010| i. The 
determination of the water deuterium fractionation in a larger 
sample of Class protostars would allow us to know, from a sta- 
tistical point of view, whether the HDO/H2O ratio is rather about 
10~ 4 - 10 3 as observed i n the NGC1333-IRAS4B proto star, in 
comets (~ 3 x 10~ 4 ; e.g. Bockelee-Morva n et al.|l998| and in 
the Earth's oceans (~ 1.5 x lQ-^ lLecuyeret al.|1998) or rather 
about a few percent, like the young stellar objects NGC1333- 
IRAS2A and IRAS 16293. With a value of about 10~ 4 in the 
protostellar phase, the fractionation ratio could be conserved 
throughout the different stages of the star formation until the 
formation of the the solar-system objects. Higher values in the 
protostellar phase invoke mechanisms in the gas phase and/or on 
the grain surfaces to explain the decrease in the deuterium frac- 
tionation for water from the protostar formation to the comets 
and solar system formation. The determination of the HDO/H2O 
ratio in a larger sample would allow us to understand whether the 
deuteration of water in protostars is somewhat similar to the so- 
lar system value and whether the environment or the conditions 
for the prestellar core phase play a role. 

Finally, we can estimate for the first time the D2O/H2O ratio 
in a low-mass protostar. According to the results of Vastel et al. 
(2010), the column density of D2O in the cold envelope (T < 
30 K) is 1.65 + 1.41 x 10 12 ctrT 2 at a 3cr uncertainty. If we con- 
sider that the absorption of the two D2O fundamental lines can 
be both due to the cold envelope (T < 30K), as well as to the ab- 
sorbing layer, the D2O/H2O ratio is on average about 10~ 3 and in 
the interval 1.1 x 10~ 4 - 3.75 x 10" 3 at 3cr. As for the D 2 0/HDO 
ratio, it is about 6% and between 0.8% and 11.6% at 3<x. The 
HDO/H 2 and D 2 0/H 2 ratio follow the statistical distri bution 
(even considering the 3cr uncertainty) as determined in |Butner] 
letaE!™' " 



(5) 



'D 2 0\ 1 /HDO\ 



These estimations of deuteration ratios in a low-mass protostar 
will certainly provide better constraints on the gas-grain chem- 
istry models. Indeed, HDO is trapped on the grain surfaces be- 
fore the thermal desorption due to the heating from the accreting 
protostar or before its photodesorption at the edges of the molec- 
ular cloud. 



5. Conclusion 

This study is the first one to use such a large number of lines 
to model deuterated water. Thanks to the numerous HDO transi- 
tions observed with Herschel/HIFI and four other lines observed 



with ground-based telescopes, we have succeeded in accurately 
determining the abundance of HDO throughout the envelope of 
the protostar and particularly the outer abundance that |Parise 
|et al.| ( |2005| l could not constrain. To estimate the abundances, 
new collisional coefficients computed with ortho and para-H2 
by Faure et al. ( 201 1 1 and for a wide range of temperatures were 
used. The best-fit inner abundance X; n (HDO) is about 1 .7 x 10~ 7 , 
whereas the best-fit outer abundance X out (HDO) is about 8 x 
10 _1 1 . To model the deep HDO absorption lines, it has been nec- 
essary to add an outer absorbing layer with an HDO column den- 
sity of 2.3 x 10 13 cirT 2 . In addition, detections of several H^O 
transitions, as well as small upper limits on other transitions, 
have allowed us to constrain both the inner and outer water abun- 
dances. Assuming a standard isotopic ratio H^O/H^O = 500, 
the inner abundance of water is about 5 x 10~ 6 and the outer 
abundance about 1.5 x 10~ 8 . The water column density in the 
added absorbing layer is about 5 x 10 14 cirT 2 . If we consider 
that this absorbing layer is created by the photodesorption of 
the ices at an Ay of ~l-4 mag, the water abundance is about 
1.5 - 3 x 10 nicely consistent with the values predicted by 
|Hollenbach et al.| (|2009). The deuterium fractionation of water 
is therefore about 1.4-5.8% in the hot corino, 0.2-2.2% in the 
colder envelope, and 4.8% in the added absorbing layer. The 3cr 
uncertainties determined in both the inner and the outer part of 
the envelope are small due to both the well constrained HDO and 
Hj 8 abundances (see Table 2). These results do not permit any 
conclusion on an enhancement of the fractionation ratio in the 
inner envelope with respect to the outer envelope. The similar ra- 
tios derived in the hot corino and in the absorbing layer suggest 
that water forms before the gravitational collapse of the proto- 
star, unlike formaldehyde and methanol, which form later after 
the CO molecules have depleted on the grains. In the cold enve- 
lope (T < 30K), the D2O/HDO ratio is estimated with a value 
of 6%. The HDO/H2O ratios found here are clearly higher than 
those observed in comets (~ 0.02%). The water deuterium frac- 
tionation ratio has to be estimated in more low-mass protostars 
to determine if IRAS 16293 is an exception among the Class 
sources or a typical protostar that could lead to the formation of 
a planetary system similar to our Solar System. In the latter case, 
processes should be at work to reprocess the water deuteration 
ratio before the cometary formation stage. 

The strong constraints obtained here emphasize the necessity 
to observe several lines in a broad frequency and energy range, 
as done with the HIFI spectral survey, to precisely estimate the 
water deuterium fractionation. In particular, the HDO lij-Oo.o 
fundamental line at 894 GHz is a key line in the modeling be- 
cause it shows emission both from the hot corino and the outer 
part of the envelope. In addition, both this transition and the lo,i- 
Oo,o transition at 465 GHz present a deep absorption probing the 
absorbing layer. The HDO 3i 2-22,i and 2 11 -2 1> 2 transitions ob- 
served at 226 and 242 GHz, respectively, with the IRAM-30m 
telescope also bring useful information as they entirely probe the 



11 



Coutens A. et aL: A study of deuterated water in the low-mass protostar IRAS 16293-2422 



hot corino. These transitions as well as several H^O transitions 
lie in the ALMA^] spectral range that will hopefully allow the 
HDO/H2O ratio to be studied and constrained in many low-mass 
protostars with very high spatial resolution. 
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